Lattice dynamics of CuGeC>3 : inelastic neutron scattering and model calculations 
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The lattice dynamics in CuGeC>3 has been analyzed by the combination of inelastic neutron 
scattering studies and lattice dynamical model calculations. We report an almost complete set 
of dispersion curves along the three orthorhombic directions and along [101]. The dispersion of 
branches associated with the modes directly involved in the spin-Peierls transition allows to explain 
the particular propagation vector of the structural distortion in the dimerized phase. 

PACS numbers: 75.40.Gb, 61.12.Ex, 63.20.-e 
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I. INTRODUCTION 

The spin-Peierls transition in CuGeOs has attracted 
considerable interest due to the relatively simple crys- 
tal structure of this material allowing experimental 
studies so £&r not possible in the organic spin-Peierls 
compoundsEI. In particular neutron scattering and 
diffraction-tpchniques could be used for the study of this 
material!. 

The spin-Peierls transition is based on magneto-elastic 
coupling : the dimerization of the magnetic chain results 
from the structural distortion. Since the gain in mag- 
netic energy is linear, and since the loss in elastic energy 
is quadratic in the distortion, in a simple one-dimensional 
picture the combined structural and magnetic transition 
must occur. In real systems inter-chain coupling may, 
however, favor three-dimensional antiferromagnetic or- 
dering. 

In CuGe03, the microscopic coupling has been clar- 
ified by experimental - magnetic and structural - and 
theoretical studies in great detail. The spin- 1/2 chains 
in CuGeOs are formed by the Cu02-chains character- 
ized by edge sharing CUO4; the oxygen in these chains 
is labeled 020. These chains are directly connected in 
the 6-direction by Ge04-tetrahedra. In contrast, the cou- 
pling along a through the apical oxygen, 01, is only weak 
due to the long CuOl-distance. A detailed description of 
the crystal structure, and its temperature dependence is 
given in references™! The variation of the magnetic in- 
teraction parameter, J, in the spin-Peierls phase results 
mainly from the modulation of the Cu-O-Cu bond-angle, 
77, which in contrast to the high,[T c superconductor par- 
ent compounds is close to 90° ml In an isolated Cu02- 
chain with 90° bond angle, there should be no antiferro- 
magnetic exchange, since two exchange-paths cancel each 
other. In CuGeOs the exact cancellation of these rather 
large values is destroyed by the deviation of the bond- 
angle from 90° and by the hybridization with the Ge- 
atoms acting as side-groups. This renders the magnetic 
interaction quite sensitive to small structural changes. 
Already in the non-dimerized phase the magneto-elastic 
coupling causes anomalous temperature dependences via 



the equilibrium of magnetic and structural energyuiJ. In 
the spin-Peierls phase the dimerization is achieved by the 
modulation of both the Cu-O-Cu bond angle 77, and the 

Sgle between the Cu02-ribbons and the Ge-Ol-bond, 
. The quantitative relation between the bond angles 
and J has been calculated by several techniques : on the 
base of band structur«-calculations0a't3, from thermody- 
namic considerationstlLand within a RPA-theory of the 
spin-Peierls transitiortj. A comparative discussion of the 
distinct techniques is given in ref.E3. 

It is obvious that phonons play an important role in 
any spin-Peierls transition. Therefore, the lattice dy- 
namics of CuGeOs has been studied by several tech- 
niques. Due to its transparency CuGeOs presents fa- 
vorable conditions for optical methods; almost all zone- 
center frequencies have been determined by the combi- 
nation of Raman and Infrared studiesE3. However, the 
spin-Peierls transition is characterized by a breaking of 
transnational symmetry, therefore, the phonons directly 
involved in the transition must have propagation vectors 
away from T. In CuGeOs the lattice is doubled along 
a and c, the propagation vector is, hence, (0.5 0.5). 
Several groups have challenged the question whether the 
involved phonon modes soften close to .the transition or 
not by inelastic neutron scatteringoOEj. These stud- 
ies tried to solve the problem by measuring only a few 
branches, which does not allow to characterize the lattice 
dynamics of this complex system. CuGe03 has 10 atoms 
in its primitive cell, the associated lattice dynamics with 
its 30 phonon branches is close to the maximum com- 
plexity one may treat today. The identification of the 
involved phonon modes could be achieved only by the 
complete study reported here. 

We have already published results concerning the 
phonon modes directly related to the spin-Peierls 
transitionO. There are four modes with the symmetry 
of the transition from Pbmm to Bbcm. Two aspects of 
these modes were unexpected. First, the distortion of the 
dimerized phase does not correspond to just one of these 
modes, but at least two polarization patterns have to be 
combined. Second, the two strongly involved modes do 
not soften close to the transition in contradiction with the 
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Cross-Fisher predictionlll. These results have stimulated 
new theoretical studies : Gros and Werner have extended 
the Cross-Fisher theory in order.-tp explain the missing 
of phonon softeningp] and UhrigcJ has stressed the im- 
portance of the high phonon frequencies compared to the 
magnetic energy scale, i.e. a non-adiabatic condition. In 
general it has been shown that coupling of the spin sys- 
tem with a branch of phonons is important for a qu.anjbi- 
tative understanding of the spin-Peierls transitionEJ'Ei 

In this paper we want to present the general lattice dy- 
namical study furnishing the basis for the identification 
of the relevant phonon modes. Only the combination 
of neutron scattering with lattice dynamical calculations 
permits to treat a phonon problem of this complexity. 
CuGe03 is one of the few complex systems where the 
lattice dynamics may be considered as being understood 
in detail. In this sense it is a rather promising material 
for validating ab-initio procedures to calculate phonon 
frequencies, so far restricted to less complex systems. In 
addition, the lattice dynamics of CuGe0 3 presents sev- 
eral aspects not related to the spin-Peierls scenario but 
interesting in themselves. Furthermore, the phonon dis- 
persion of CuGeOa should have some exemplary charac- 
ter for the wide class of silicates and germanatescJ and 
also for other cuprate chain systems of current interest. 

In section II we describe the neutron scattering stud- 
ies and present the lattice dynamics of CuGe03. The 
phenomenological model used to accompany the mea- 
surements and to interpret the data is introduced in sec- 
tion III. The general discussion of the phonon dispersion 
in CuGeOa based on the experimental studies and the 
model is given in section IV. In section V we give some 
further information on the modes directly involved in the 
transition which is not included in our previous paper. 
Finally, in section VI, we briefly discuss the possibility of 
distinct structural order parameters. 



II. INELASTIC NEUTRON SCATTERING 
STUDIES 

The main part of the studies of the phonon disper- 
sion has been performed on the triple axis spectrometer 
ITa, operated by the Forschungszentrum Karlsruhe at 
the Orphee reactor in Saclay. A few additional studies 
have been made on_the spectrometer 2T also installed at 
the Orphee reactorcf . 

Since it has been necessary to study the phonon dis- 
persion in CuGe03 in its entity, the amount of beam time 
spent on this problem was considerable. 25 days on IT 
and 7 days on 2T were used in order to determine an 
almost complete set of phonon branches along the three 
orthorhombic directions and along [101]. For the study 
of the temperature dependence of modes relevant to the 
spin-Peierls transition, another 8 days were used on ITa. 

For a first experiment, two crystals of 550 mm 3 volume 
each have been coaligned in the [100]/[001]-orientation 
using a double goniometer. Due to the long 6-axis this 



orientation allows attaining Q-values out of the scattering 
plane by tilting the sample with the goniometer. Along 
the [100]-direction almost all modes could be determined 
in this orientation. One of the two crystals shows a 
larger mosaic spread around a, which was tolerable in 
the [100]/[001]-orientation since it becomes hidden by 
the large vertical divergence. For orientations with a not 
parallel to the scattering plane it has been advantageous 
to use only one crystal. Further experiments were per- 
formed in the [010]/[001]-, [100]/[010]- and [101]/[010]- 
orientations. 

For all measurements a pyrolytic graphite (PG) ana- 
lyzer was used; in the low frequency range up to ^lOTHz 
a PG-(002) monochromator and, at high frequencies, 
a Cu-(lll) monochromator. For a few studies requir- 
ing an exceptional resolution a Cu-(220) monochroma- 
tor was mounted. Due to double focusing arrangement 
of monochromators and analyzers the_spectrometer IT 
yields a considerable gain in intensity^. The focusing 
arrangement is incompatible with the use of collimators 
but yields a better resolution in comparison to a configu- 
ration with flat crystals and open collimations. Roughly, 
the focusing configuration corresponds to a conventional 
spectrometer with 37'-collimators throughout. All scans 
were performed with the final energy fixed in the neutron 
energy loss mode and, apart a few scans aiming at acous- 
tic phonons, the final energy was fixed to the values of 
3.555 and 7.37THz, where higher order contaminations 
may be suppressed by the PG filter. 

The phonon dispersion was studied mainly at room 
temperature; just for the highest frequency modes it has 
been favorable to cool the sample to 10K in order to 
reduce the background. 

In figure 1 we show typical phonon scans which illus- 
trate the different effort to be made to analyze phonon 
frequencies at low and high energy. The intensity, J, of a 
onp-phonon process in neutron energy loss mode is given 
byH 



/ oc !.(n(a,)+l){V -^ = . e (-^Q- d ).(Q. ed )} 2 (1), 



{u> denotes the phonon frequency, q the wave-vector, 
Q = g + q the scattering vector (g the reciprocal lattice 
vector) e~ Wd the Debye- Waller-factor, and the sum ex- 
tends over the atoms in the primitive cell with mass 
scattering length bd, position r<i and polarization vector 
ed. In general, the intensity is determined by the Bose- 
factor, n(ui), and the ^-term, which strongly reduce the 
effectiveness to observe high energy modes. In analogy 
to the elastic structure factor, the sum in (1) is called dy- 
namic structure factor; it describes the interference of the 
interaction of the distinct atoms weighted by the scalar 
product of their displacements with the scattering vec- 
tor Q. Only phonons with some polarization parallel to 
Q can be observed. The right side of equation (1) may 
be calculated with a lattice dynamical model in order to 
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14 15 16 17 18 19 22 24 26 
v (THz) v (THz) 

FIG. 1: Exemplary scans aiming at determining phonon fre- 
quencies at the T-point and in the Brillouin-zone with fits by 
Gaussian profiles and sloping background. 



III. LATTICE DYNAMICAL MODEL 
CALCULATIONS FOR CUGEO3 

In order to analyze a complex problem like the phonon 
dispersion it is essential to take profit of symmetry con- 
siderations. The 10 atoms in the primitive cell cor- 
respond to 30 zone-center frequencies characterized by 
their polarization vectors. The symmetry of the CuGe03 
crystal structure allows a-separation corresponding to ir- 
reducible representationsEH each of them being charac- 
terized by a certain polarization scheme. The T-modes 
may be divided into : \A g + 2A U + B lg + 6Bi u + 4B 2g + 
6£?2u + 3-B3 9 + 4i?3 U . The schemes of the polarization 
patterns are given in Table I; they result from the crys- 
tal symmetry. A g modes show the full symmetry, their 
displacement parameters hence correspond to the free pa- 
rameters in the CuGeC>3 structure, Ge-x, Ol-x, 02-x and 
02-y. A u -modes show displacements of Cu and 02 along 
c; they are silent since shifts in neighboring chains cancel 
each other. Bi U are characterized by displacements of all 
atoms in the i-direction with in-phase shift of equivalent 
atoms; in particular they contain the starting points of 
the acoustic branches. i?i S -modes are even modes which 
break a two- fold axis. 



predict favorable conditions for the measurement and in 
order to identify single phonon modes. 

Figure 1 shows a scan across a transverse acoustic 
phonon at g=(0 2) (upper left); the mode is polarized 
along the c-direction. This mode yields a high peak inten- 
sity in accordance with the strong (0 2) Bragg-reflection 
intensity. In contrast, the determination of the highest 
frequencies requires an extreme effort in beam-time; the 
figure 1 (lower right) presents results for a strong dy- 
namic structure factor; modes with unfavorable structure 
factors may not be analyzed in practice in this frequency 
region. In the medium energy range the analysis of the 
phonon spectra is dominated by the question whether the 
modes may be isolated for a certain Q-value. 

The phonon frequencies were obtained by fitting Gaus- 
sian distributions to the measured profiles. This proce- 
dure may be used in principle only when the curvature 
of the dispersion surface can be neglected in the range of 
the four-dimensional resolution ellipsoid. The whole sot 
of dispersion curves along the orthorhombic directions, 
given in figure 2, is separated according to the irreducible 
representations. Only a few high energy branches were 
not studied in our experiments; however, their zone cen- 
ter frequencies are fixed by the optical techniques and 
the dispersion is most likely flat. 



4 




0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 

(x 0) (x 0) (x 0) (x 0) (0 x 0) (0 x) (0 Ox) (0 x) (0 x) 



FIG. 2: Phonon dispersion in CuGeOs: branches are separated according to irreducible representations along the three or- 
thorhombic directions; circles denote experimental points (open circles along (0x0) Si- and E2-modes) and lines the frequencies 
calculated by the model. 
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TABLE I: Polarization schemes according to the crystal structure of CuGeOs for all r-modes and for the representations along 
the three orthorhombic directions and along [101] and for the Tj-modes associated with the spin-Peierls transition at q=(0.5 
0.5); the first lines give the positions of 10 atoms forming a primitive unit, the following lines show the displacements of these 
atoms. A letter at the i-position, signifies that this atom is moving along the i-direction, a second appearance of the same 
letter signifies that the second atom moves with the same amplitude ("-" denotes a phase shift) in the corresponding direction. 
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10 Ax : 4 A 9 + 6 B lu 10 A 3 : 4 B 2g + 6 B 2u 

3 A 2 : 2 A u + Bi g 7 A 4 : 3 B 39 + 4 B 3tl 

10 Si : 4 A 9 + 6 B 2 „ 10 £ 3 : 4 B 2 g + 6 B lu 

5 S 2 : 2 A„ + 3 B 3g 5 E 4 : B lg + 4 B 3tl 

8 Ai : 4 A 9 + 6 B 3u 9 A 3 : 3 B 3 g + 6 Bi„ 

6 A 2 : 2 A u + 4 B 2g 7 A 4 : Big + 6 B 2tl 

17 XOXj : 8 Aj + 9 A 3 : 4 A 9 + 6 B 3u + 3 B 3a + 6 Bi. 

13 X0X 2 : 6 A 2 + 7 A 4 : 2 A„ + 4 B 2ff + Bi fl + 6 B 2u 



TABLE II: Compatibility relations in CuGc0 3 

In an analogous way the crystal symmetry allows also 
to divide the phonon polarization patterns for any q- 
value at high symmetry points or lines in the Brillouin- 
zone, in particular along the three orthorhombic direc- 
tions: A corresponds to [xOO], £ corresponds to [0x0] 
and A corresponds to [00x]. The separation of the 30 
branches according to four distinct representations is an 
important simplification for the analysis of the phonon 
dispersion as it may already be seen in figure 2. 

In the definition used here, the Xi-representation 
always contains the longitudinal acoustic modes with 
X=A, £ or A. The ^-representation is characterized 
by the fact that no acoustic modes correspond to it, and 
X%- and ^-representations are given by the transverse 
acoustic branches in the sequence of the polarizations, 
along a before along b before along c. The polarization 
patterns according to these representations are given in 
Table I for the main symmetry directions and for [X0X]. 
The division according to the representations has signifi- 
cance for the neutron scattering intensity : they yield se- 
lection rules for the observation of the modes which allow 
an identification. The relations between the irreducible 
representation at T and in the zone are given by the 
compatibility relations shown in Table III. For instance 
Ag-modes always correspond to the Xi-representation 
and the A u to the ^-representation. Since CuGeC>3 
is an insulator, polar modes exhibit Lydane-Sachs- Teller 
(LST) splitting, the longitudinal frequencies correspond 
to starting points of Ai-branches with X being the po- 
larization direction, and the transverse frequencies to the 
X3 and Xj-branches. The symmetry further leads to a 
degeneration at q=(0 0.5 0), where a Ei-branch connects 
with a S3-branch and where a X^-branch connects with 
a E4-branch. 

In the frame of harmonic lattice dynamics one may 
reduce the equations of movement to a 3n-dimcnsional 
Eigen-valu£.._problem for each allowed q-value, see 
referencesEElo, w 2 e = De; here e is the 3n-dimensional 
Polarization vector and D the dynamical matrix given by 



D a , (d, d!) = J /2 ^(Od, l'd')exp(iql') (2). 



Here the indices I' numerate the cells, d,d' the atoms 
within one cell, and a, (3 the three space directions. 
^a,/3(0d, I'd!) are the force constants between the atoms 
d and d! in cells shifted by V corresponding to the direc- 
tions a and (3. The determination of the force constants 
is the central problem to analyze the phonon dispersion. 

For the description of the phonon dispersion in 
CuGeC>3 we use Coulomb-potentials with effective 
charges, V(r) oc Zl z r 2e and the repulsive forces are de- 
scribed by Born-Mayer potentials V(r) = A-exp(—r/r ). 
Shell charges describing a single-ion polarizability have 
been introduced for all atoms with a single shell-core 
force constant. The interatomic, forces act on the shells 
in our model. Chaplot et alHd have reported a com- 
mon model for compounds related to the high-T c super- 
conducting cuprates; these Cu-O-parameters have been 
used as starting values for the description of CuGeOs. 
The model has been initially adapted to the optical fre- 
quencies and has been continuously refined with the in- 
elastic neutron scattering results. The parameters have 
been fitted to the observed phonon frequencies and to 
zero forces on the atoms in the equilibrium positions. In 
addition, the predictions of the model concerning the dy- 
namic structure factors have been compared to the mea- 
sured intensities. AlL-ealculations were performed with 
the GENAX progranM 

We found that several features of the phonon disper- 
sion could only be reproduced by the inclusion of angular 
forces for the Ge-O-interaction. The angular forces, 4-r, 
are divided by the lengths of the two distances, in order 
to yield values comparable to usual force constants. The 
angular forces at the Ge-site reflect the strong covalent 
character of these bonds and are expected. However, also 
an angular force in the Cu-O-arrangement gave a signif- 
icant improvement, whereas, in the HTSC cuprates no 
such parameters are needed. 

The angular forces, however, interfere with the forces 
arising from the Born-Mayer potentials. Therefore, it is 
impossible to describe all identical pairs by a single po- 
tential. In a first model, model I see Table III, additional 
force constants were introduced for nearest neighbors, the 
forces acting in this model correspond hence to the sum 
of the direct force constants plus the contribution from 
the Born-Mayer potentials. The Cu-02 potential as well 
as the van der Waals terms were completely dropped in 
this model. Model I gives a very good description of the 
frequency data : the entire set of over 700 phonon fre- 
quencies is described with a mean deviation of 0.19 THz. 
We then attempted to simplify the model by limiting the 
amount of force constants, model II in Table III. Only 
nearest neighbor Ge-0 interaction is described by force 
constants, all other forces are deduced from potentials. In 
addition 0-0 van der Waals terms and another angular 
force were introduced. This model too gives a satisfac- 
tory description of the frequency data, mean deviation 
of 0.21 THz, but the the equilibrium forces on the atoms 
are much higher. 
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ionic part 
model I model II 



ion 


Z 


Y 


K 


Z 


Y 


K 


Cu 


1.77 


4.0 


2.8 


1.87 


3.91 


2.0 


Ge 


2.26 


0.0 


/ 


1.836 


0.0 


/ 


01 


-1.22 


-2.8 


2.0 


-1.03 


-2.18 


1.8 


02 


-1.41 


-3.3 


2.0 


-1.34 


-3.14 


1.8 





potentials 








model I 


model II 


pair 


A(eV) r (A) 


A(eV) 


ro (A) 


Cu-Ol 


900 0.305 


1008 


0.305 


Cu-02 


/ / 


4294 


0.226 


Ge-0 


2500 0.243 


/ 


/ 


Ol-Ol 


1300 0.288 


2000 


0.284 


01-02 


1800 0.288 


2000 


0.284 


02-02 


1500 0.288 


2000 


0.284 


force constants (dyn/ 


:m) 




pair 


F G 


F 


G 


Ge-02 


98846 30082 


603041 


-53026 


Ge-Ol 


2233 31133 


470090 


-37047 


Cu-02 


262013 -29617 






02-02 (in Cu0 4 ) 


-15963 4619 






01-02 


5202 -4591 






02-02 (in Ge0 4 ) 


16428 -8400 






angular force constants (c 


yn/cm) 




Ol-Ge-02 


2523 


8776 




Cu-02-Ge 


1705 


4936 




Ge-Ol-Ge 


10816 


19754 




02-Cu-02 


3753 


3921 




02-Ge-02 


/ 


10897 





TABLE III: Model parameters for the description of the 
phonon dispersion in CuGeOs; for the explanation of the 
parameters see text; Z,Y are in electron charges; K in 
10 6 dyn/cm; in model II the uniform O-O potentials contains 



also a van der Waals term -30(ev A 6 
atomic distance in A. 



with r the inter- 



Both CuGe03-models present ionic charges much 
lower than the chemical values, as it is typically ob- 
served. The strong reduction, at least in case of Ge, is 
certainly caused by strong covalence of the bonds. The 
shell charges and forces, however, could be chosen similar 
to values found in other materialsEJ. Due to its better 
description, we discuss only the results of model I in the 
following. 





Popovic et al. 


neutrons 


Udagawa et al. 


calculation 




5.606 


5.532 


5.516 


5.81 




9.953 


9.692 


9.923 


9.92 




17.808 


18.13 


17.78 


17.85 




25.753 


/ 


25.723 


25.69 


Big 


11.632 


11.68 


11.512 


11.68 


B 2 g 


3.478 


3.46 


3.388 


3.29 




6.716 


6.86 


6.626 


6.62 




12.801 


12.95 


12.92 


12.85 




26.352 


26.71 


26.382 


25.80 


B 3g 


3.298 


3.441 


3.298 


3.65 




12.322 


12.34 


12.292 


12.82 




21.345 


/ 


21.376 


22.14 



TABMC IV: Comparison of thei Raman-results by Udagawa 
et alE3 and by Popovic et al.E3 with the frequencies deter- 
mined by inelastic neutron scattering and calculations with 
the lattice dynamical model, all values are given in THz. 



IV. DISCUSSION OF THE PHONON 
DISPERSION IN CUGEO3 



groups have performed Raman-studies on 
he agreement between these results is 

4 



Several 
CuGeQ^j 

goodlHH&K-i r-Iii Table IV we compare the results of 
Udagawa et al.EHl and Popovic et alJla with the values 
obtained by inelastic neutron scattering and the frequen- 
cies calculated with the model. The experimental values 
agree within 2% and also the agreement with the model is 
satisfactory, except for the Ge-O-bond stretching vibra- 
tions. The latter failure is certainly due to the insuffi- 
cient description of the covalent Ge-O-bonds; within the 
parameter range a better description may be obtained 
but at expense of the low frequency agreement. 

Infra*cistjudji^r«pCuGe03 were reported by several 
groupsll§BEal23'EI(HL3, whose results correspond well for 
modes polarized along b and c. The determination of the 
frequencies for modes polarized along a is hampered by 
the shape of the crystal and was attempted only by one 
groupEl Table V compares the optical results with the 
neutron and model frequencies; there is only one discrep- 
ancy concerning an a-polarized mode. 

Popovic et al. have interpreted a weak signal near 
480cm~ 1 or 14 THz as a i?i u -mode. This frequency is, 
however, not in accordance with the lattice dynamical 
calculations. The determination of the LST-splitting by 
inelastic neutron scattering is in general difficult, since 
the extension of the four-dimensional resolution ellipsoid 
always favors the observation of the transverse modes. 
For this reason there is little chance to determine the LO- 
frequency if the LST-splitting is comparable or smaller 
than the resolution. Therefore, we have attempted the 
measurement of the longitudinal polar frequencies only 
for a few favorable cases. Figure 3 shows an energy scan 
across the B3„ pair at 15.7-17.3 THz; these frequencies 
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Popovic et al. TO/LO 


neutrons 


calculation 




3.927 4.017 


4.16 / 


4.081 4.212 




9.833 11.153 


9.3 9.6 


9.252 9.312 




14.33 14.66 


9.47 10.91 


9.860 11.294 




18.438 19.182 


18.24 / 


17.945 19.052 




24.07 25.27 


// 


24.960 25.800 


B2u 


1.439 1.469 


1.508 / 


1.370 1.438 




6.326 6.925 


6.372 / 


6.335 7.023 




8.544 9.234 


8.58 / 


8.961 9.419 




11.302 12.412 


11.41 / 


11.03 12.156 




23.145 25.693 


// 


23.42 25.068 


B3„ 


4.947 5.037 


5.109 / 


5.239 5.249 




15.869 17.990 


15.73 17.27 


15.435 17.77 




21.586 24.074 


// 


22.15 23.23 




/ 


3.081 


2.811 




/ 


15.42 


14.77 



TABLE V: Comparison of the infrared-results by Popovic et 
al.E3 with the inelastic neutron results and calculated frequen- 
cies, TO and LO denote the transverse and longitudinal optic 
frequencies, respectively, in THz. 
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FIG. 3: Scans across polar modes: above B3„ pair, below 
Biu-pair. 

confirm the infrared results. In the lower scan pee ex- 
pects the Bin-mode, assumed by Popovic et alJlj near 
14THz. Unambiguously, this scan reveals a lower energy 
for this mode, in good agreement with the model calcu- 
lations. From this scan and from the extrapolation of 
the branches starting at this mode, we conclude that the 
TO- and LO-frequencies of the third highest -Bi„-mode 
amount to 9.47 and 10.91 THz. The next lower B Xu - 
mode is found at 9.3-9.6 THz, hence there is an overlap 
of these two pairs, which will produce a single plateau 
in the reflectivity spectra at 9-11 THz, in addition this 



mode possesses little oscillator strength. Indeed we. may 
describe the spectrum observed by Popovic et al£3 rea- 
sonably well on the basis of these frequencies. 

The crystal symmetry only yields the schemes of the 
polarization patterns as given in Table I. For any repre- 
sentation with a multiplicity higher than one, one has to 
know the force constants in order to determine the po- 
larization patterns. Figure 6 shows the polarization pat- 
terns calculated with the lattice dynamical model. Sim- 
ilar results were also reported by Popovic et afij; small 
deviations are probably due to the misinterpretation of 
the -Biu-mode. 

The A^-modes correspond to displacements of the Cu 
and 02 positions along c. The mode with an opposite 
shift causes significant alternation of the Cu-02-bond dis- 
tance, this mode has a high frequency of 15.4 THz. The 
in-phase vibration, where the entire ribbons are shifted 
along c, is low in frequency at 3.1 THz. The distortion 
in the spin-Peierls phase is related to both j4„-modes. 

The rotation of the Cu02-ribbons around the c-axis 
corresponds to the j4 9 -mode at 9.9 THz. This high fre- 
quency may surprise in view of the fact that the ribbons 
rotate as function of temperature around d3. The polar- 
ization pattern of the 9.9 THz-mode lacks the alternation 
of the 6-lattice constant and the associated shift of the 
Ge04-tetrahedra; without these elements the vibration 
is hence rather hard. 

Several modes correspond to rotations of the tetrahe- 
dron chains around the c-axis. The i?2g-niode at 6.7 THz 
may be described by the rotation around an axis near the 
middle of the 02-02-edges; this vibration is rather hard. 
The -B2g-mode at 3.5 THz, in contrast, corresponds to the 
rotation around the Ol-Ol-line and is hence related to 
the distortion in the spin-Peierls phase. In this low fre- 
quency i?2g-mode neighboring tetrahedra rotate in the 
same sense, and the corresponding anti-phase rotation 
is approximatively realized in the -B2u-niodes at 6.3 and 
1.5 THz. The pure anti-phase tetrahedron rotation is 
not an Eigen-mode due to the stronger influence of the 
CuO-forces; the difference between the two i?2u modes 
concerns just the Cu-shift. 

Modes with dominant Cu-contribution have lower fre- 
quencies due to the higher mass; frequencies are partic- 
ularly low for displacements perpendicular to the Cu-O- 
bonds. In the lowest optic B\ u mode, Cu atoms neighbor- 
ing along b shift almost in the direction of their connec- 
tion. The mode with the opposite phase corresponds to 
the lowest -B2u-niode with a significantly lower frequency, 
1.4 THz. This pattern allows the 02-atoms to follow the 
Cu-atoms yielding an almost rigid shift of the ribbons. 
Also, the tetrahedra are only slightly distorted in this 
pattern, since they rotate around the c-axis. In accor- 
dance rwith the two-dimensional character of the CuGeOa 
latticeQ this vibration may be considered as the trans- 
verse vibration of the zick-zack-planes. This explains the 
exceptional low frequency of this mode which is related to 
the low lying longitudinal acoustic Si branch, see below. 

The compared experimental and calculated phonon 
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FIG. 4: Polarization patterns of all F-modes and their fre- 
quencies in THz. The figure on the top illustrates the crystal 
structure of CuGeC>3 and gives the polarization pattern of 
the zone boundary mode of the longitudinal acoustic branch 
(q=(0 0.5 0) ). 



dispersion shown in figure 2 indicates strong dispersion 
mostly in the low frequency range. Furthermore, the 
curves are steeper along c than along & or a. This re- 
flects the arrangement of the strong bonds, since the 
phonon dispersion may be considered as the Fourier- 
transformation of the force constants. The layered struc- 
ture of CuGeC>3 has no covalent bonds along a, therefore, 
the dispersion stays flat in this direction, at least in the 
higher frequency range. In consequence, one may expect 
that also the character of the vibrations does not change 
in this direction. In contrast, the chain configuration 
along c will mix the character in the Brillouin-zone along 

[ooi]. 

The dispersion along b is determined by the coupling 
between CuOvribbons and tetrahedron chains. The 
glide mirror plane in space group Pbmm transforms into 
the second formula unit in the primitive cell; for vanish- 
ing coupling one might separate the zone center modes 
into modes with the two units vibrating in- or out-of- 
phase but the finite coupling results in a dispersion along 
b. The degeneration of two branches at (0 0.5 0) corre- 
sponds in most cases to the connection of such a pair. 
The flat dispersion of many of the branches along b shows 
that the coupling is weak for those modes. 

In the following we discuss now several particularities 
of the phonon dispersion in CuGeC>3 which are not di- 
rectly related to the spin-Peierls transition. 

- Flat longitudinal acoustic branch along b and associ- 
ated modes - Lorenzo et al. have reported a longitudinal 
acoustic branch in the 6-directiontj, which has been in- 
terpreted as being essential for the spin-Peierls transition. 
The low-lying,b*anch has been qualitatively confirmed in 
later studiesEjO and also by our own results. However, 
the results do not agree quantitatively. The low-lying 
branch may be described within the lattice dynamical 
model. The low frequencies are explained due to the 
lowest B2u-niode. From this mode a Si branch starts, 
which interacts with the longitudinal acoustic branch of 
the same symmetry. The character of the acoustic vi- 
bration is hence transferred to the optic branch, and the 
acoustic branch carries the character of the B2 M -mode at 
1.5THz, as it is illustrated by the polarization pattern of 
the zone-boundary acoustic mode, see top of figure 4. 

The flat LA branch is hence the consequence of the low- 
lying optical mode. Similar effects are also seen in the 
other directions where the lowest B2„-mode belongs to 
the A3- or A4-representations. The corresponding branch 
interacts with the transverse acoustic branches, A3 and 
A4, respectively. 

The flattening was interpreted asjiidicating direct rele=. 
vance for the spin-Peierls transitions. Damascelli et al.o 
report for the lowest B2 U -mode an anomalous softening 
upon cooling, which we do not confirm. The polariza- 
tion pattern of these vibrations, however, has no simi- 
larity with the spin-Peierls distortion in CuGeC>3. The 
low frequencies just reflect a structural instability of this 
material. 

The agreement between frequencies obtained by differ- 



10 




0.25 0.50 0.5 1.0 1.5 

(0x0) v (THz) 

FIG. 5: a) Frequencies of the longitudinal acoustic mode in 
the S-direction; the measured frequencies of this work were 
corrected for resolution-effects. For comparison we also show 
results of other groupst-J'Ej. The dashed line corresponds to 
the linear extrapolation from the elastic constant determined 
by Poirier et al.cil, and the solid line to the frequencies calcu- 
lated with our model, b) Scan calculated by the convolution 
of the dispersion surface with the spectrometer resolution for 
Q=(0 6.15 0); the exact value of the LA mode frequency at 
this q- value is indicated by the broken line. 
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ent groups is poor for the LA-Ei-branch, see figure 5. 
The convolution of the experimental resolution with the 
dispersion surface may explain these differences. Figure 
5b) shows a scan calculated by the convolution of the ex- 
perimental resolution (including the mosaic spread of the 
sample) with the calculated dispersion. By several cal- 
culations of this type the resolution-induced correction 
factors have been determined and applied to the mea- 
sured frequencies. The corrected frequencies are shown 
in figure 5. Resolution effects are most likely the origin of 
the bad agreement between the different neutron results. 

- Soft-mode-behavior - Figure 6 shows the low en- 
ergy part of the phonon dispersion of the A4- and XOX2- 
branches. Due to the larger amount of branches of the 
same symmetry, the interpretation is rendered difficult 
in the [101]-direction, but one may recognize that the 
branch starting at the B2„-mode at 6.3 THz softens con- 
tinuously through the zone, exhibits several interactions, 
and ends finally at the lowest zone-boundary frequency. 
The a-component does not seem to be relevant for this 
decrease. As shown in figure 4, the B2„ mode corre- 
sponds to a displacement of 02-positions in a square in 
the same direction coupled with a movement of the Cu- 
atoms. Along the c-direction the pattern changes contin- 
uously, for q=(0 0.5), 02-02-edges neighboring along c 
are moving in opposite directions, this pattern is shown 
in figure 7, it may be characterized as a folding of the 
ribbons. In the following this mode is called soft mode 
(SM). 

The frequency decrease in these branches can be con- 
sidered as soft-mode behavior; it points to some struc- 
tural instability, for instance this mode softens by 2.3% 
upon temperature decrease from 295 to 12 K. This insta- 
bility in CuGeC>3 is still far from condensing, but upon 
substitution or pressure it might become enhanced. The 



FIG. 6: Part of the dispersion of the A4- and X0X2-branches. 



SM-instability, however, is not directly related to the 
spin-Peierls transition. 

The A4-mode of similar frequency corresponds to a Cu- 
shift, see figure 7. The mode with opposite phase between 
neighboring chains is the A3-mode, in contrast to the 
corresponding T-modes, Bi„ and ^>2u, the chain coupling 
only weakly affects the zone-boundary modes. 

- Behavior of the highest Ge-O-branches - The disper- 
sion of the highest Ge-O-branches is not fully understood. 
In general the description of these branches by our model 
is poor due to the insufficient description of the covalent 
forces, see above. The behavior of the highest A2- and 
XOA^-branches, shown in figures 2 and 11, is surpris- 
ing. Again, the influence of the a-component is negligi- 
ble. The bending down of the branch towards the zone- 
boundary is reproduced by the model only in general. 
We do not think that this behavior results from magneto- 
elastic coupling, but it appears more likely to be an in- 
trinsic property of the tetrahedron chains, "ab-initio" 
calculations should give more insight into this problem. 

- Further predictions of the lattice dynamical model - 
Table VI compares the elastic constants determined by 
ultra sound and by Brillouin scattering techniques with 
those calculated by the model. The first ultrasound jnea- 
surements were limited to the diagonal termscirCJO and 
a complete set of elastic constants was detepnined only 
by Ecollivet et al. using Brillouin scatteringCj. There is 
some scattering in the experimental data obtained with 
the ultra-sound techniques. The model *alues compare 
best to the combined study in referencaU. The lattice 
model predicts C22 to be significantly higher than the 
experimental value, as may be also seen in figure 6. The 
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FIG. 7: Polarization-pattern of the four lowest modes at q=(0 
0.5). 

ell c22 c33 c44 c55 c66 cl2 cl3 c23 
US-a 0.624 0.372 3.264 ////// 
US-b 0.74 0.21 3.32 ////// 
US-c 0.66 0.345 2.79 ////// 
US-d 0.71 0.345 3.43 0.37 0.33 0.22 III 
BR 0.64 0.376 3.173 0.353 0.353 0.184 0.321 0.469 0.227 
model 0.823 0.500 3.457 0.408 0.435 0.165 0.297 0.403 0.223 

TABLE VI: Comparison of the experimentally determined 
elastic constants (ultra-sound : US-a Peirier et al.Lj, US-b 
Poirier et al£3, US-c Saint-Paul et al.,Ej, US-idn Ecollivet et 
al.J-j; Brillouin scattering : BR Ecollivet et al.,Lj) with those 
calculated with the lattice dynamical model; elastic constants 
are given in 10 12 dyn/cm 2 . 



model seems to slightly underestimate the interaction be- 
tween the optic and acoustic branches near the zone cen- 
ter. Upon increase of the weight of the LA modes one 
may obtain better agreement for the elastic constants 
but at the expense of the description at higher energies, 
which has been considered to be more relevant. 

The phonon density of states, PDOS, lias been stud- 
ied by Arai et alJlj and by Fujita et al.c£l by inelastic 
neutron scattering on polycrystalline samples. In figure 
8 we show the PDOS weighted with the atomic cross 
section divided by the mass, the generalized PDOS (GP- 
DOS), calculated with the lattice dynamical model; the 
agreement with the experiment is at most qualitative. 
However, the good agreement of the model with the op- 
tic frequencies and with the dispersion curves determined 
here excludes the experimental GPDOS. We think that 
the deviation results from insufficient correction of multi- 
phonon processes. Also the ..strong temperature depen- 
dences point to such effectsE§EJ. Neither our studies, nor 
the-pptical experiments nor the similar work by Nishi et 
al.tj give evidence for a strong temperature dependence 
of the phonon spectrum in CuGeOa. In addition, the 
measured GPDOS presents huge intensity above the cut- 



FIG. 8: GPDOS calculated with the lattice dynamical model. 



off value of 26 THz, which however, is already fixed by 
the optics studies. 

From the calculated PDOS one may obtain the spe- 
cific heat at constant volume. Most of the published 
work on the specific he a t i de als with the anomaly at,the 
spin-Peierls transitionE3c30'ED. Only Weiden et al.c3 re- 
port a measurement up to room temperature which is 
compared in figure 9 with the results of the calculations. 
The difference between specific heat at constant volume 
and at constant pressure has been corrected empirically, 
indicated by the solid and dashed lines in figure 9. Since 
the coefficient of thermal expansion enters the difference 
in second order, the correction is extremely small, it re- 
mains almost in the width of the lines. Surprisingly, the 
measured specific heat is always about 12% higher than 
the calculated one; at least at higher temperature this 
can not be attributed to magnetic contributions. Again, 
it seems unlikely that the lattice dynamical model is that 
wrong; for instance the large number of high optic fre- 
quencies unambiguously indicate that the specific heat at 
room temperature is far from its saturation value corre- 
sponding to the Dulong-Petit law. A larger specific heat 
might arise from strong anharmonicity, as-.in the classic 
example of Agl where a sub-lattice meltsEa. The anhar- 
monicities in CuGeOs appear to be, however, too small 
for such an explanation. 

Another check of the model may be obtained through 
the comparison with the experimental Debye- Waller fac- 
tors, see Table VII. The agreement at room temperature 
is very good, and also the low temperature coefficients, 
which are more difficult to determine, arc well described!! 



V. MODES IN RELATION TO THE 
SPIN-PEIERLS TRANSITION 

The structural part of the spin-Peierls transition is 
characterized by a symmetry reduction in the crystal 
structure from Pbmm to Bbcm. The relevant phonon 
modes may be identified with the hela of group theory. 
In the notation by Stokes and HatcfH these modes are 
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FIG. 9: Comparisonrejf the experimental specific heat (points, 
taken from reference!— I) with the calculated ones; the solid line 
corresponds to c p and the dashed one to cy . 



atom exp. cal. 


atom exp. cal. 


Ge-Ull 70(3) 55 
-U22 101(3) 90 
-U33 40(3) 31 

Cu-Ull 115(3) 100 
-U22 142(2) 116 
-U33 46(3) 42 
-U12 50(2) 37 


01- Ull 81(4) 64 
-U22 129(4) 117 
-U33 48(3) 42 

02- U11 144(3) 127 
-U22 178(3) 161 
-U33 61(3) 56 
-U12 79(2) 69 



TABLE VII: Debye-Waller-coefficients in CuGe0 3 at room 
temperature : compariscm of experimental (average of x-ray 
and neutron results fromn) and calculated values, in 10 _4 A 2 . 



labeled T+ at q=(0.5 0.5), TheT+-modes are charac- 
terized by four independent parameters in their polariza- 
tion pattern which reflect the four additional structural 
parameters in the dimerized structure. The correspond- 
ing irreducible representation is one-dimensional and has 
a multiplicity of four. The T ^-modes have been discussed 
in detail in referencealaO including their temperature 
dependence. In this section we want to give some fur- 
ther information on the dispersion of the branches con- 
nected to these four modes directly involved. The lowest 
T^-mode at 3.12 THz is B2 ff -like associated with a rota- 
tion of the 02-02-edges around the c-axis; it modulates 
the 02-02-Ge angle, S, and, therefore, also the magnetic 
interaction J. The second lowest mode at 6.53 THz is 
associated with the modulation of the Cu-02-Cu bond 
angle r/, which determines strongest the magnetic inter- 
action. This 77-mode seems to be the one most important 
for the spin-Peierls transition. The higher frequency T^"- 
modes at 11 and 25 THz correspond to the Cu-02- and 
Ge-0 bond stretching and are less relevant for the spin- 
Peierls transition. The T^~-representation corresponds to 
the propagation vector of q=(0.5 0.5), therefore, it is in- 
teresting to analyze the dispersion in the [101]-direction, 
see figure 11. 

- [101]- direction and identification of the relevant 
modes - The T-modes whose polarization patterns are 
closest to the structural distortion in the spin-Peierls 
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FIG. 10: Polarization patterns of the two modes with Cu- 
shift along c at q=(0 0.5); the upper two pictures show a 
projection of a single Cu04-square. 
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FIG. 11: Dispersion of the branches in the [101]-direction. 
The symbols denote measured frequencies and the solid lines 
those calculated by the lattice dynamical model. Dashed lines 
denote the dispersion of corresponding A-branches neglecting 
the influence of the a*-component. 



phase are the A u and B2 9 -modes where branches of A2- 
symmetry start in the A-direction. As may be seen 
in the polarization schemes given in Table I, the A2- 
representation allows the z-shift of Cu and the x, y-shifts 
for 02. These Cu and O-characters of these modes are 
well separated near the zone center but they mix in the 
zone. Along the [101]-direction the modes may be divided 
only according to two different representations, 17 X0X\- 
and 13 X 0X2-modes. The large number of modes of the 
same symmetry renders the interpretation of the data 
rather difficult. The compatibility relations show that 
the T^-modes correspond to the zone-boundary modes 
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of the AOAVbranches, see Table I. 

At the zone-boundary (0 0.5) there are two vibrations 
where the Cu moves along c, the upper one is isolated in 
frequency near 12 THz and the lower one near 7.7 THz. 
As may be seen in figure 10, the difference in these modes 
consists in the phases of Cu and 02-displacements. In the 
higher mode, Cu moves towards the shorter 02-02-edge 
yielding a strong alternation of the Cu-O-bond lengths 
and, therefore, the higher energy. In contrast, for the 
lower mode Cu shifts towards the longer edge yielding an 
alternation of the bond angle. Besides the a-component 
this movement represents the main part of the structural 
distortion in the dimerized phase. These modes are called 
77-modes in the following. 

The B2 9 -mode near 3.5 THz corresponds to the rota- 
tion of the edges of the CuC-4-squares and hence to the 
second dominant feature of the spin-Peierls distortion. 
A flat A2-branch is starting at this mode, which ends 
at the zone-boundary at a mode with opposite twisting 
of edges neighboring along c. Again, only the a compo- 
nent is missing for the description of the distortion in the 
dimerized phase. In general one may arrive at the rele- 
vant T^-modes by either analyzing the path from (0 
0.5) to (0.5 0.5) where there is only little dispersion.ar 
studying the [101]-dispersion, see below and referenceO. 

Figure 11 shows the phonon dispersion in the [101]- 
direction. Again the dispersion is well described by the 
model. From the compatibility relation in Table II one 
recognizes that AOXi-branches correspond to the combi- 
nation of Ai- and A3-branches in the [001]-direction, and 
AOAVbranches correspond to A2- and A4-branches. In 
figure 11 we show the corresponding A-branches as bro- 
ken lines neglecting the x-component of the wave-vector 
(i.e. in a projection). The resemblance shows that the 
a*-component influences the frequencies only slightly, in 
agreement with the flat dispersion along a and the lack 
of any covalent bond in this direction. However, the in- 
fluence of the a*-component may not at all be neglected 
for the interpretation of the intensities. Along the A- 
direction one observes only Ai-branches at Q=(0 2+x); 
the intensities of A3-modes are suppressed, similar to an 
elastic extinction rule. In the [101]-direction, however, 
one may also observe the AOAi-modes corresponding to 
the A3-modes at Q=(x 2+x) with comparable inten- 
sities. The observation that structure factors are much 
more sensitive than the intensities corresponds to the ex- 
perience with small structural distortions^. 

Due to the larger number of modes of same symme- 
try the mixing of the characters is even more important 
in the [101]-direction. Nevertheless, it has been possible 
to clcarlvJdentify the T^-modes at bfq=(0.5 0.5) in 
referenceO. The T^-mode rotating the 02-edges is found 
at 3.12 THz, the most relevant ^-modulating mode at 
6.53 THz, and the two highest T^-modes exhibit frequen- 
cies of 11.2 and 24.9 Thz, respectively, see referenceO. 

- Comparison with Raman-studies in the dimerized 
phase- Raman-studies below Tsp Jaave been very 
early performed by Udagawa et al.E3 and by Sugai 



et al. U: they show additional peaks in the spin- 
Peierls phase. In the meanwhile there have been 
many publications on these intensities, which agree 
concerning t he .. experim e ntal, Ji ndings but not on the 
interpi^atiore^BME9k3ka>H&a^. All groups find 
five additional peaks; additional much weaker intensities 
are reported in reference! 5 and appear to arise from 
some different mechanism. 

The lowest intensity near 30cm is of magnetic ori- 
gin. All Raman-interpretations agree that the highest 
two strong peaks correspond to folded phonons, they 
perfectly agree to our results on the two highest Tg - 
modes. But also the two peaks at lower frequencies which 
were interpreted controversely may be identified as T^"- 
modes, since their frequencies perfectly agree with those 
of the B2 S -like and 77- T^~-modes. The neutron scatter- 
ing studies reveal a strong broadening of the 77-mode in 
good-agreement with the Raman-experiments; Lemmens 
et al£9 report a width of 16cm _1 . In addition, it has been 
found that just the Raman-intensity corresponding to the 
77-mode can be observed a few degrees K above the spin- 
Peierls transition l 6C U 61 l This underlines once more the 
exceptional role of the ?7-modc in the spin-Peierls transi- 
tion. 

- Dispersion of the branches connected to the - 
modes - In general the dispersion is rather flat along the 
a-direction. This is also valid for branches starting at 
the T^-modes and passing along the path (0.5-x 0.5) 
on the zone boundary. The endpoint corresponds to the 
zone-boundary of the A-direction. The small influence of 
the a* component is already illustrated in figure 11 and 
was used to discuss the character of the T^-modes. The 
modes at q=(0 0.5) may be further divided according 
to symmetry . In particular, there is a representation 
of multiplicity four which corresponds to the T^-modes 
at (0.5 0.5). The frequencies of these four modes lie 
at 3.29, 7.69, 12.1 and 24.9 THz, compared to the T£- 
frequencies at 3.12, 6.53, 11.2 and 24.9 THz at room tem- 
perature. 

The (0 0.5)-mode corresponding to the B2 9 -like Tj- 
mode turning the 02-edges is only slightly harder than 
the T^-mode, at room temperature. At low temperature, 
the difference is even smaller since the (0 0.5)-mode 
hardens slightly to 3.34 THz, whereas the T^-mode is sig- 
nificantly shifted due to the magneto-elastic couplingO. 

The difference in frequency for the modes modulating 
the Cu-O-Cu-bond angle, 77, the main feature of the spin- 
Peierls distortion in CuGeOa , is more pronounced, about 
1 THz. This dispersion explains the doubling of the lat- 
tice in the a-direction; the distortion with modulation 
along a just requires less structural energy. For these 77- 
modes too, the stronger frequency increase upon cooling 
is observed for q=(0.5 0.5); but also the (0 0.5)-mode 
frequency increases significantly, by 1.5% . 

The second highest T^-mode exhibits a significant fre- 
quency increase towards (0 0.5), whereas the Ge-O- 
modes are little influenced by the a*-component. 

The dispersion along the 6-direction on the zone- 
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FIG. 13: Polarization patterns ol four modes at q=(0.5 0.5), 
which are related to the T^-modes at 6.6 THz (77-mode) and 
at 3.1 THz (B 29 -like mode). 



FIG. 12: Dispersion of branches connected to the T^-modes 
along the path (0.5 x 0.5) at room temperature ; symbols 
denote experimental frequencies and lines calculations. 



boundary (0.5 x 0.5) is comparable to that in the S- 
direction. In particular, there is a degeneration of two 
modes at (0.5 0.5 0.5). Again, two (0.5 0.5)-modes con- 
nect which correspond to the phase shift between the two 
formula units. Figure 12 shows the dispersion along this 
path, and figure 13 indicates the polarization patterns of 
several (0.5 0.5)-modes relevant for the following dis- 
cussion. 

There are four modes in the range 3-4 THz at q=(0.5 
0.5): the T2 -mode, two Cu-modes, associated to the 
lowest optic Bi u - and B2 U -frequencies and the soft mode. 
As in the E-direction, the two Cu-modes connect. For the 
T^-mode the situation is more complicated : consider- 
ing the scheme of a single Cu02-ribbon, one might expect 
an opposite-phase mode, where along b neighboring 02- 
02-cdges rotate in the same sense concerning their angle 
to the a-axis. This movement, however, would result in 
a strong modulation of the Ge-O-bonds causing much 
higher frequency. On the contrary one may consider this 
T+-mode as a rotation of the Ge04-tetrahedra around 
the Ol-Ol-line. In the T^-mode two tetrahedra neigh- 
boring along b turn both clockwise or both anti-clockwise; 
the movement with different phase - one clockwise the 
other anti-clockwise - results in a completely different 
displacement scheme concerning the Cu02-ribbons, since 
both 02-positions in an 02-02-edge perpendicular to c 
move then in the same direction. This pattern corre- 
sponds to the soft mode at (0.5 0.5), see figure 13. The 
lowest T^-mode connects hence with the soft mode along 
the 6-direction but looses the modulation of the magnetic 
interaction along this, path. The anomalous frequency in- 
crease upon coolingEZI observed for the T^"-mode has been 
followed along b*; the frequency hardening continuously 
decreases till (0.5 0.5 0.5); in contrast the soft-mode even 



softens by 2.3(3) % upon cooling from room temperature 
to 12 K. 

For the 77-mode the Ge-O-forces do not play an impor- 
tant role; therefore, one may easily follow these modes 
along b* . One finds a significant frequency enhancement; 
the corresponding mode is found at 8.56 THz, its po- 
larization pattern, see figure 13, exhibits a smaller Cu- 
displacement than that of the T^^-mode. In the range 
5-6 THz, there are several flat branches included in figure 
12 which are not related to the spin-Peierls distortion. 
The temperature dependence of the bond angle modu- 
lating modes close to (0.5 0.5) could not be studied, 
since in this region additional magnetic intensity appears 
upon coolingE-3. For qi,=0.5 no temperature dependence 
has been detected, also the intensities near 5.8 THz are 
almost temperature independent. 

The dispersion of the branches connected to the two 
lower T^-modes explains the propagation vector of the 
structural distortion in the spin-Peierls phase, q=(0.5 
0.5). The 6-component is fixed through the Ge-O-forces, 
which shift the mode rotating the 02-02-edges to very 
high energies and which also case a sizeable slope of the rj- 
modulating branch. The weaker inter-plane forces lead to 
the doubling of the a-parameter mainly due to their influ- 
ence on the 77-modes. The structural distortion involved 
in the spin-Peierls transition has, therefore, a clear three- 
dimensional character. Therefore, a RPA treatment of 
the spiiwahonon coupling in the transition appears to be 
justified!! 

The dispersion of the phonon branches should reflect 
the extension of the critical scattering along the three or- 
thorhombic directions. In a classical soft-mode transition 
one may expand the frequency of the involved mode at 
the propagation vector of the structural distortion, qo, 
here (0.5 0.5) with q = q + q a + qb + q c : 



w(q) = uj + c a • q a 2 + c b • q b 2 + c c • q c 2 (4). 
Linear terms do not exist, since there must be a fre- 
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flections which are not sensitive to the same features of 
the distortion. Intensity at (0.5 3 0.5) is mainly related 
to the distortion corresponding to the lowest T^-mode 
with B2 g -character, and intensity at (0.5 4 0.5) reflects 
the ry-modulation. In spite of their different sensitivities 
these two super-structure exhibit exactly the same tem- 
perature dependence as it is shown in figure 14. Fitting a 
critical behavior, / cx (Tsp—T) 213 , to the temperature de- 
pendence of these superstructure intensities yields a criti- 
cal exponent .ofpff= 0.30(1) in good agreement with previ- 
ous resultsE3iSI 'E3 and in particular with the thermal ex- 
pansion studyo. The scaling between the two superstruc- 
ture reflections remains also valid slightly above the spin- 
Peierls transition, indicating that the whole structural 
transition may be described by a single order-parameter. 



FIG. 14: Intensity of the two superstructure reflections at 
Q=(0.5 3 0.5) and (0.5 4 0.5) as function of temperature. The 
two intensities reflect approximately the distortion of the two 
lowest T^-modes. The line corresponds to a fit by power law. 

quency minimum at qo. The ratio of the constants 
c a ,Cb,c c should correspond to the ratio of the correla- 
tion lengths of critical scattering in the three directions. 
In CuGeC>3 this analysis is not obvious since there is no 
softening and since at least the correlation length along c 
is dominated by the magnetic interaction. Furthermore, 
the fitting of the experimental dispersion by a quadratic 
expansion is not very satisfactory; one obtains c a = 2.6 
THz/A 2 , c h =8.0 THz/A 2 and c Q /c 6 = 0.3, for the 77-mode 
which is the most relevant one. 

The observed correlation lengths decrease very rapidly 
above the transition; only a few degrees K above it they 
amount to a few lattice constants; therefore, rather large 
q a and qb-values are involved, where a quadratic expan- 
sion is certainly no longer valid. Nevertheless, the ratio 
between the correlation lengths determined by Schocffcl 
et al.E3 |^ ~ 1/4 may be qualitatively explained by the 

dispersion of the branches around the T^-77-mode. 

VI. POSSIBILITY OF INDEPENDENT ORDER 
PARAMETERS 

Since the structural distortion in the spin-Peierls phase 
does not correspond to a single phonon Eigen-mode, one 
might speculate whether the low temperature phase is 
described by a single order parameter or whether more 
than one parameters are needed in order to describe the 
temperature dependence of the distortion. We have ana- 
lyzed this problem by measuring two superstructure re- 



VII. CONCLUSION 

The various inelastic neutron scattering studies have 
yielded a detailed understanding of the lattice dynamics 
in CuGeC>3. The developed model not only describes the 
observed dispersion curves, but it allows the interpreta- 
tion of Raman and infrared-frequencies, elastic constants 
and anisotropic Debye- Waller parameters. It further pre- 
dicts the PDOS and the phononic part of the specific 
heat, on which only little experimental information ex- 
ists till today. 

Only the extensive study of the phonon dispersion has 
allowed the identification of those modes directly involved 
in the transitiontll. The dispersion of the branches con- 
necting with these modes clearly explains the occurrence 
of the structural transition at q=(0.5 0.5), since the 
phonon frequencies of modes modulating the bond-angle 
77 are lowest at this q-value. The pronounced dispersion 
of the phonon branches involved in the spin-Peierls tran- 
sition clearly illustrate a three-dimensional character of 
the structural part of the transition. The influence of the 
magneto-elastic coupling on the temperature dependence 
of phonon frequencies should occur not only in CuGeOa 
but also in materials with related structure, like for ex- 
ample the spin-ladder compounds. 

CuGeOa is amongst the most complex materials where 
the phonon dispersion has been analyzed in that detail. 
It appears therefore to be a promising candidate for the 
extension of current "ab initio" techniques to complex 
materials. 
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